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Abstract. We investigate the star formation history and 
chemical evolution of low surface brightness (LSB) disk 
galaxies by modelling their observed spectro-photometric 
and chemical properties using a galactic chemical and 
photometric evolution model incorporating a detailed metal- 
licity dependent set of stellar input data. For a large frac- 
tion of the LSB galaxies in our sample, observed properties 
are best explained by models incorporating an exponen- 
tially decreasing global star formation rate (SFR) ending 
at a present-day gas fraction M &as / (M gas + M sta rs) = 0.5 
for a galaxy age of 14 Gyr. For some galaxies small am- 
plitude star formation bursts are required to explain the 
contribution of the young (5-50 Myr old) stellar popu- 
lation to the galaxy integrated luminosity. This suggests 
that star formation has proceeded in a stochastic manner. 

The presence of an old stellar population in many late- 
type LSB galaxies suggests that LSB galaxies roughly fol- 
low the same evolutionary history as HSB galaxies, except 
at a much lower rate. In particular, our results imply that 
LSB galaxies do not form late, nor have a delayed onset 
of star formation, but simply evolve slowly. 

Key words: Galaxies: Low Surface Brightness spirals - 
Galaxies: formation and evolution - Galaxies: fundamen- 
tal parameters 



1. Introduction 

Deep searches for field galaxies in the local universe have 
revealed the existence of a large number of galaxies with 
low surface brightnesses that are hard to detect against the 
night sky (see the reviews by Impey & Bothun 1997 and 
Bothun, Impey & McGaugh 1997). Many different kinds of 
LSB galaxies have so far been found, including giant LSB 
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galaxies (e.g. Sprayberry et al. 1995) and red LSB galax- 
ies (O'Neil et al. 1997), but the most common type seems 
to be "blue LSB galaxies" : late-type, disk-dominated spi- 
rals with central surface brightnesses Ho{B) ^ 23 mag 
arcsec~ 2 . 

In this paper, we will concentrate on these late-type 
LSB spirals. Observations show that they are neither ex- 
clusively dwarf systems nor just the faded counterparts of 
"normal" high surface brightness (HSB) spirals (de Blok et 
al. 1996, hereafter dB96). In many cases, LSB galaxies fol- 
low the trends in galaxy properties found along the Hubble 
sequence towards very late types. These trends include in- 
creasingly blue colours (e.g. Ronnback 1993; McGaugh & 
Bothun 1994; de Blok et al. 1995 [hereafter dB95]; Bell et 
al. 1999), decreasing oxygen abundances in the gas (e.g. 
McGaugh 1994; Ronnback & Bergvall 1995; de Blok k 
van der Hulst 1998a), and decreasing Hi surface densities 
(from type Sc onwards; e.g. van der Hulst et al. 1993, here- 
after vdH93). Despite the low gas densities, LSB spirals 
rank among the most gas-rich disk galaxies at a given to- 
tal luminosity as their Hi disks in general are extended 
(Zwaan et al. 1995; dB96). The fact that LSB galaxies 
still have large reservoirs of gas together with their low 
abundances suggest that the amount of star formation in 
the past cannot have been very large (vdH93; Van Zee et 
al. 1997). Clearly, LSB galaxies are not the faded remnants 
of HSB spirals. 

The unevolved nature of LSB spirals can be interpreted 
in various ways. For instance, LSB spirals could be systems 
whose stellar population is young and in which the main 
phase of star formation is still to occur. Alternatively, the 
stellar population in LSB spirals could be similar to those 
in HSB galaxies but with a young population dominating 
the luminosity. This has been explored and modelled in 
studies by e.g. Knezek (1993), Jimenez et al. (1998), dB96, 
Gerritsen & de Blok (1999) and more recently in two pa- 
pers by Bell et al. (1999) and Bell & de Jong (1999). These 
studies all indicate a scenario in which LSB galaxies arc 
unevolved systems with low surface densities, low metal- 
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licities and a constant or increasing star formation rate 
(cf. Padoan ct al. 1997, however, for a different opinion). 

In this paper we address these various scenarios for the 
evolution of LSB galaxies by detailed modelling of their 
spectro-photometric and chemical properties. Model re- 
sults are compared directly with the observed colours, gas 
phase abundances, gas contents, and current star forma- 
tion rates of LSB galaxies. Using simple model star for- 
marion histories we confirm the results of previous studies 
and show that the ratio of young stars to old stars in LSB 
galaxies is larger than typically found in HSB galaxies. 

In reality, however, star formation is unlikely to pro- 
ceed as smoothly as the various population synthesis mod- 
els suggest. Gerritsen & de Blok (1999) point out that 
small surges in the star formation rate are likely to be 
very important in determining the observed colours of LSB 
galaxies. In this paper we will pay special attention to the 
influence of small bursts of star formation on the optical 
colours, and show that they can have a significant influ- 
ence on the inferred properties of LSB galaxies. 

This paper is organized as follows. We briefly compare 
observational data of LSB galaxies with those of spirals 
and dwarf galaxies in Sect. 2. In Sect. 3, we describe the 
ingredients of the galactic evolution model we used. Some 
general properties of the model are presented in Sect. 4. 
The comparison with observational data is presented and 
discussed in Sect. 5. The impact of small amplitude star 
formation bursts in LSB galaxies is investigated in Sect. 
6 and predicted star formation rates are compared with 
the observations in Sect. 7. A short discussion of where 
LSB galaxies fit in the grand scheme of galaxy formation 
and evolution is presented in Sect. 8. A short summary is 
given in Sect. 9. 

2. Observational characteristics of LSB galaxies 

2.1. Sample selection 

We use the sample described in dB96. This consists of 24 
late- type LSB galaxies (inclinations up to ~ 60°), taken 
from the lists by Schombert et al. (1992) and the UGC 
(Nilson 1973). The sample is representative for the LSB 
galaxies generally found in the field. We selected a sub- 
sample of 16 LSB galaxies for which good data are avail- 
able. For these systems, optical data have been taken from 
dB95, Hi data from dB96, and abundance data from Mc- 
Gaugh & Bothun (1994) and de Blok & van der Hulst 
(1998a). 

The galaxy identification and Johnson U BV and Kron- 
Cousins RI are given in columns (1) to (6) in Table 1 (a 
Hubble constant of Ho = 100 km s _1 Mpc -1 was used). 
Column (7) lists the -B-band luminosities; column (8) the 
Hi-mass-to-light ratio; column (9) and (10) the Hi and dy- 
namical masses respectively. The gas fraction fi rot in col- 
umn (11) is given by M gas /(M gas + M* jlnax ) where M* iInax 



denotes the mass of the stellar component obtained from 
maximum disk fitting of the rotation curve (dB96) (see 
Sect. 2.6). The mean [O/H] abundances are listed in col- 
umn (12). 



2.2. Magnitudes and colours 

We compare in Fig. la and b magnitudes and broadband 
colours of LSB galaxies to those of late-type HSB face-on 
spirals (de Jong & van der Kruit 1994) and dwarf galaxies 
(Melisse & Israel 1994; Gallagher & Hunter 1986, 1987). 
The usual distinction between HSB and LSB galaxies is 
purely an artificial one since normal galaxies along the 
Hubble sequence show a continuous range in central sur- 
face brightness, i.e. ranging from values around the Free- 
man value for early type systems to the very faint values 
observed for the late- type LSB galaxies in our sample (e.g. 
de Jong 1995; dB95). For the purpose of this paper we 
will continue to use the distinction as a shorthand, where 
"HSB galaxy" should be interpreted as "typical Hubble se- 
quence Sc galaxy." Fig. la and b show that LSB galaxies 
are generally bluer and fainter than their HSB counter- 
parts (e.g. McGaugh 1992; vdH93; dB95). 



2.3. Abundances 

Estimates of the ISM abundances in LSB galaxies pre- 
dominantly rely on abundance determinations of their Hn 
regions. Oxygen abundances of Hn regions listed in Ta- 
ble 1 are taken from McGaugh (1994) and de Blok & van 
der Hulst (1998a). 

We assume that the Hn-region abundances on average 
are a reasonable indicator of the ISM abundances within 
a given LSB spiral. The intrinsic scatter in [O /H] among 
different Hn regions within a given LSB galaxy is usually 
of the order 0.2 dex (e.g. McGaugh 1994). Bell & de Jong 
(1999) find that the gas metal abundances in LSB galaxies 
generally trace the stellar metallicities quite well, with the 
stars being ~ 0.4 dex more metal-poor. 

In Fig. lc we compare mean [O /H] abundances of Hn 
regions in LSB galaxies with those in HSB spirals (Zarit- 
sky et al. 1994) and dwarf galaxies (Melisse & Israel 1994; 
Gallagher & Hunter 1986, 1987). On average, LSB galaxies 
generally follow the correlation between the characteristic 
gas-phase abundance and luminosity as found for HSB spi- 
rals (Zaritsky et al. 1994), but with a large scatter. The 
range in [O /H] at a given B magnitude is nearly one dex. 
This scatter is probably related to evolutionary differences 
among the LSB galaxies of a given B magnitude (e.g. in 
the ratio of old to young stellar populations) and/or the 
Hn regions they contain. 
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Fig. 1. Observational data on LSB galaxies compared to that for normal spirals and dwarf galaxies. Symbols refer to 
LSB galaxies (triangles; data from de Blok et al. 1995, 1996), normal face-on spirals (full dots; de Jong & van der 
Kruit 1994), and dwarf galaxies (crosses; Mellise & Israel 1994). LSB galaxies which are probably dwarf systems are 
indicated by triangles with crosses overlayed. Typical error bars are indicated in the bottom left of each panel. 
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Table 1. Observational data on LSB galaxies (de Blok et al. 1995, 1996) 
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— 18.0 


— 18.4 
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2.0 
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F568-3 
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F571-5 


-16.6 


-16.5 


-16.9 


-17.1 


* 


8.79 


* 


9.99f 


10.45 


* 


-1.53 


F571-V1 


* 


-16.4 


-16.9 


-17.3 


* 


8.75 


1.2 


8.82 


10.15 


0.33 


-0.94 


F574-2 


* 


-17.0 


-17.7 


-17.8 


* 


8.99 


0.9 


8.97 


9.51 


0.66 


* 


F577-V1 


-17.9 


-17.6 


-18.0 


-18.1 


* 


9.23 


0.8 


9.15 


9.20t 


* 


* 


U0128 


-18.5 


-18.2 


-18.7 


-18.9 


-19.3 


9.48 


1.2 


9.56 


10.86 


0.39 


* 


U0628 


-18.5 


-18.5 


-19.1 


-19.3 


-19.7 


9.59 


* 


* 


10.69 


* 


* 


U1230 


-18.9 


-17.7 


-18.2 


-18.5 


-18.9 


9.28 


1.7 


9.50 


10.80 


0.60 


-0.76 


Typ. Dwarf 


-17.5 


-17 


-17.5 


-18 


-18.5 


9.0 


1.0 


9.0 


9.0 


0.7 


-0.6±0.4 


Typ. LSBG 


-18 


-17.5 


-18 


-18 


-18.5 


9.2 


1.3 


9.3 


10.0 


0.5 


-0.6±0.4 


Typ. HSBG 


-20 


-19.5 


-20 


-20.5 


-21 


10.0 


0.4 


9.6 


11.0 


0.1 


-0.0 



* not available, f uncertain 



2.4- Extinction 

Observational studies by (amongst others) Bosnia et al. 
(1992), Byun (1992), Tully & Verheijen (1997) show that 
late-type spirals appear transparent throughout their disks. 
This supports the idea that face-on extinction in LSB 
galaxies is relatively low, i.e. Eb-v Si 0.1 (e.g. McGaugh 
1994). This is consistent with findings by Tully & Verhei- 
jen (1998), who studied a large sample of galaxies in the 
UMa cluster and found the LSB galaxies in their sample 
to be transparent, whereas non-negligible extinction was 
present in the HSB galaxies. 

In this study we will ignore the possible effects of dust 
on the colours and magnitudes. See Sect. 5.3 for a discus- 
sion. 



2.5. Gas masses 

We compare in Fig. Id the observed atomic gas masses 
M g ~ 1.4Mhi (corrected for helium) of LSB galaxies with 
those of HSB galaxies and dwarfs. Strictly speaking we 
should also take the molecular gas mass into account. 
However, so far no CO-emission has been detected in late- 
type LSB galaxies (Schombert et al. 1990; de Blok & van 
der Hulst 1998b, but see Knezek 1993 for a discussion 
on giant LSB galaxies). The amount of molecular gas in 
LSB galaxies is probably small, even though relatively 
high CO/H 2 conversion factors may apply in these low 
metallicity galaxies (Wilson 1995, Mihos et al. 1999). We 
conclude that, on average, LSB galaxies are considerably 



more gas-rich (up to a factor — 3) than HSB galaxies of 
the same B magnitude (see also Fig. lg). 

2.6. Gas fractions and mass-to-light ratios 

We show in Figs, lg and h the distribution of the mass-to- 
light ratio Mhi/ Lb vs (B — V) and vs fi rot , respectively. 
LSB galaxies exhibit considerably higher Mm/ Lb ratios 
than HSB spirals (cf. Table 1 and McGaugh & de Blok 
1997). 

We would however like to know the ratio of gas mass 
to stellar mass, i.e. the gas fraction fix = M gas / (M gas + 
Mstars of a galaxy. This involves the conversion of the 
observed galaxy luminosity to stellar mass. However, the 
mass-to-light ratio of the underlying stellar population is 
generally not well known and, in fact, is an important 
quantity to determine. Assuming a fixed value of the mass- 
to-light ratio will however introduce artificial trends since 
this ratio is expected to vary among galaxies having differ- 
ent star formation histories (e.g. Larson & Tinslcy 1978). 

As a partial solution we have determined the gas frac- 
tion ^ rot = M gas /(M gas + M* jmax ) where M» iinax denotes 
the mass of the stellar component obtained from maxi- 
mum disk fitting of the rotation curve (dB96). This method 
likely overestimates the contribution of the luminous stel- 
lar disk to the observed total mass distribution by about a 
factor of 2 (e.g. Kuijken & Gilmore 1989; Bottema 1995) 
and, therefore, also provides a lower limit to the actual 
gas fraction. In summary, while \i\ is a model quantity, 
^ ro t is determined from observations. 
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Figs, le and If show the distribution of the present- 
day gas fraction [i rot and total gas mass vs (B — V). It can 
be seen that LSB galaxies and dwarfs exhibit larger gas 
fractions than HSB spirals. Taking into account the factor 
of ~ 2 overestimate, values for LSB galaxies are typically 
/Ui ^ 0.5 and ~ 0.1 for HSB spirals. 

We conclude that the available evidence supports of 
the view that LSB galaxies are relatively unevolved sys- 
tems compared to HSB spirals. This agrees well with the 
fact that LSB galaxies usually display properties interme- 
diate to those of HSB spirals and dwarf galaxies. 

3. Model description and assumptions 

We here summarise the galactic evolution model used (for 
a more extensive description see van den Hock 1997, here- 
after vdH97). We concentrate on the stellar contribution 
to the total galaxy luminosity in a given passband (other 
contributions are neglected). For a given star formation 
history (SFH), we compute the chemical enrichment of a 
model galaxy by successive generations of evolving stars. 
To derive the stellar luminosity in a given passband at 
a given age, we use a metallicity dependent set of theo- 
retical stellar isochrones as well as a library of spectro- 
photometric data. 



as input for the spectro-photometric evolution model de- 
scribed below. 

Boundary conditions to the chemical evolution model 
are the galaxy total mass M to t, its evolution time t cv , 
and the initial gas abundances. Note that solutions of the 
galactic chemical evolution equations are independent of 
the ratio of the SFR normalisation and M to t- Primordial 
helium and hydrogen abundances are adopted as Y p = 
0.232 and X = 0.768 (cf. Pagcl & Kazlauskas 1992). Initial 
abundances for elements heavier than helium are initially 
set to zero. 

We list the main input parameters in Table 2, i.e. the 
adopted IMF-slope, minimum and maximum stellar mass 
limits at birth as well as the progenitor mass ranges for 
stars ending their lives as AGB star, SNIa, and SNII, re- 
spectively. For simplicity, we assume the stellar yields of 
SNIb,c to be similar to those of SNII. Furthermore, we 
assume a fraction jy SNIa = 0.015 of all white dwarf pro- 
genitors with initial masses between ^2.5 and 8 M Q to 
end as SNIa. These and other particular choices for the 
enrichment by massive stars are based on similar models 
recently applied to the chemical evolution of the Galactic 
disk (e.g. Groenewegen et al. 1995; vdH97). 

3.2. Spectro-photometric evolution model 



3.1. Chemical evolution model 

We start from a model galaxy initially devoid of stars. We 
follow the chemical enrichment during its evolution assum- 
ing stars to be formed according to a given star formation 
rate (SFR) and initial mass function (e.g. a power law 
IMF: dN /dm = M(m) cx m 7 ). Both stellar and interstel- 
lar abundances as a function of galactic evolution time t ev 
are computed assuming that the stellar ejecta are returned 
and homogeneously mixed to the ISM at the end of the 
individual stellar lifetimes (i.e. relaxing the instantaneous 
recycling approximation; sec Searle & Sargent 1972). A 
description of the set of galactic chemical evolution equa- 
tions used can be found in e.g. Tinsley (1980) and Twarog 
(1980), see also vdH97. This model is a closed-box model. 

We follow the stellar enrichment of the star forming 
galaxy in terms of the characteristic element contributions 
of Asymptotic Giant Branch (AGB) stars, SNII and SNIa. 
This treatment is justified by the specific abundance pat- 
terns observed within the ejecta of these stellar groups 
(see e.g. Trimble 1991; Russell & Dopita 1992). A detailed 
description of the metallicity dependent stellar lifetimes, 
element yields, and remnant masses is given by vdH97 
and van den Hoek & Groenewegen (1997). We compute 
the abundances of H, He, O, and Fe, as well as the heavy 
element integrated metal-abundance Z (for elements more 
massive than helium), during the evolution of the model 
galaxy. Both the SFH, IMF and resulting element abun- 
dances as a function of galactic evolution time are used 



The total luminosity of a galaxy in a wavelength interval 
A A is determined by: 1) the stellar contribution L„, 2) 
the ISM contribution L; sm (e.g. Hn-regions, high-energy 
stellar outflow phenomena, etc.), and 3) the contribution 
due to absorption and scattering L cxt - 

^gal(^) = L* + Lism — Lcxt (1) 

where each term in general is a complex function of galac- 
tic evolution time. We concentrate on the stellar contribu- 
tion and neglect the last two terms in Eq. (1). Then, the 
luminosity in a wavelength interval AA at galactic evolu- 
tion time t = T can be written as: 

L^(t = T)= (2) 

nT rm (T-t) 

\ I L^ x (m,Z(t),T -t)S(t)M(m)dmdt 

where m\ denotes the lower stellar mass limit at birth, 
m (t) the turnoff mass for stars evolving to their remnant 
stage at evolution time t ev , and L^ x the luminosity of a 
star with initial mass m, initial metallicity Z(t), and age 
(T-t). Wc assume a separable SFR: S{m, t) = S(t)M(m) 
where S(t) is the star formation rate by number [yr _1 ] and 
M(m) the IMF [M Q -1 ]. By convention, we normalise the 
IMF as J M(m)dm = 1 where the integration is over the 
entire stellar mass range [mi, m u ] at birth (cf. Table 2). 

Starting from the chemical evolution model described 
above, we compute the star formation history S(m,t), 
gas-to-total mass-ratio fj,(t), and age-metallicity relations 
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Table 2. Parameters related to IMF and stellar enrichment 

7 —2.35 slope of power-law IMF 

(mi, m u ) (0.1, 60) Mq stellar mass range at birth 

(m* 013 , m^ GB ) (0.8, 8) M Q progenitor mass range for AGB stars 

(mf NH , mjf" 1 ) (8, 30) Mq progenitor mass range for SNII 

(mf NIa , mf" a ) (2.5, 8) M© progenitor mass range for SNIa 

i/ SNIa 0.015 fraction of progenitors ending as SNIa 



(AMR) Zi(t) for different elements i. Thus, at each galac- 
tic evolution time t ev the ages and metallicities of pre- 
viously formed stellar generations are known. To derive 
the stellar passband luminosity A we use a set of the- 
oretical stellar isochrones, as well as a library of spectro- 
photometric data. Stellar evolution tracks provide the stel- 
lar bolometric luminosity effective temperature T e g, 
and gravity g, as a function of stellar age for stars with ini- 
tial mass m born with metallicity Z». We compute Eq. (2) 
using a spectro-photometric library containing the stellar 
passband luminosities L^ x tabulated as a function of T e g, 
log 5, and Z» (see below). 

The turnoff mass m (T—t) occuring in Eq. (2) depends 
on the metallicity Z(t) of stars formed at galactic evolu- 
tion time t ev . For instance, the turnoff mass for stars born 
with metallicity Z = 10~ 3 at a galactic age of t cv = 14 
Gyr is m ~ 0.8 M Q (e.g. Schaller et al. 1992). This value 
differs considerably from m ~ 0.95 Mq for stars born 
with metallicity Z = Zq . Such differences in m affect the 
detailed spectro-photometric evolution of a galaxy by con- 
straining the mass-range of stars in a given evolutionary 
phase (e.g. horizontal branch) at a given galactic evolu- 
tion time. In the models described below, we explicitly 
take into account the dependence of m (t) on the initial 
stellar metallicity Z„ (see vdH97) 

3.3. Stellar evolution tracks and spectro-photometric data 

We use the theoretical stellar evolution tracks from the 
Geneva group (e.g. Schaller et al. 1992; Schaerer et al. 
1993). These uniform grids cover large ranges in initial 
stellar mass and metallicity, i.e. m = 0.05 — 120 Mq and 
Z = 0.04 — 0.001, respectively. These tracks imply a re- 
vised solar metallicity of Zq = 0.0188 with Y Q = 0.299, 
and AY/AZ = 3.0 for a primordial He-abundance of 
Y p = 0.232. For stars with m > 7 Mq, these tracks were 
computed until the end of central C-burning, for stars with 
m = 2 - 5 Mq up to the early- AGB, and for m < 1.7 M Q 
up to the He- flash. For stars with m ^ 0.8 Mq, we used the 
stellar isochrone program from the Geneva group (Maeder 
& Meynet, priv. comm.; see also vdH97 ). 

To cover the latest stellar evolutionary phases (i.e. hor- 
izontal branch (HB), early- AGB, and AGB) for stars with 
m ^ 8 Mq, we extended the tracks from Schaller et al. 
with those from Lattanzio (1991: HB and early-AGB) for 
m <~ 1 — 2 Mq, and from Groenewegen & de Jong (1993; 
early-AGB and AGB) for m ~ 1-8 Mq. These tracks 



roughly cover the same metallicity range as the tracks 
from the Geneva group. Corresponding isochrones were 
computed on a logarithmic grid of stellar ages, covering 
galactic evolution times up to t cv <~ 14 Gyr. Isochrones 
are linearly interpolated in m, log Z , and log t. 

The spectro-photometric data library that we use is 
based on the Revised Yale Isochrones and has been de- 
scribed extensively by Green et al. (1987). These data in- 
clude stellar UBVRI Johnson-Cousins magnitudes cov- 
ering the following ranges in T eff [K] = 2800 to 20000, 

log g [cm s~ 2 ] = -0.5 to 6, and log {Z/Zq) 2.5 to 

+0.5. Corresponding spectro-photometric data for stars 
with T c ff > 20000 K have been adopted from Kurucz 
(1979) at solar metallicity, covering T off = 20000 - 50000 
K. 

Although a detailed description of the tuning and cal- 
ibration of the adopted photometric model is beyond the 
scope of this paper, we note that the model has been 
checked against various observations including integrated 
colours and magnitudes, luminosity functions, and colour- 
magnitude diagrams of Galactic and Magellanic Cloud 
open and globular clusters covering a wide range in age 
and metallicity (see vdH97). 

4. General properties of the model 

Before discussing the modelling of LSB galaxies, we first 
describe the general behaviour of the models assuming a 
few different star formation scenarios. 

We start from a model galaxy with initial mass M g (t = 

0) = 10 10 Mq, initially metal-free and devoid of stars. 
The chemical and photometric evolution of this galaxy 
is followed during an evolution time t cv = 14 Gyr, as- 
suming one of the theoretical star formation histories dis- 
cussed below. Unless stated otherwise, we assume that 
stars are formed according to a Salpeter (1955) IMF (i.e. 
7 = —2.35) with stellar mass limits at birth between 0.1 
and 60 Mq (cf. Table 2). 

The following star formation histories are considered: 

1) constant, 2) exponentially decreasing, 3) linearly de- 
creasing, 4) exponentially increasing, and 5) linearly in- 
creasing. Normalised SFRs and resulting agc-ISM metal- 
licity relations are shown in Fig. 2. For each model, the 
amplitude of the SFR is chosen such that a present-day 
gas-to-total mass-ratio [i\ = 0.1 is achieved. 

In columns (2) to (6) of Table 3, we list the functional 
form of the SFR, average past and current SFRs, the ra- 
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Fig. 2. Theoretical star formation histories (left) and resulting [O/H] vs age relations (right): constant SFR (thick 
solid line), linear decrease (thin solid), exponential decrease (dashed; r s f r = 5 Gyr), fast exponential decrease (dotted; 
T s f r = 4 Gyr), linear increase (dash-dotted), exponential increase (dot-dashed). SFRs have been normalised to a current 
gas-to-total mass-ratio zzi =0.1. 



Table 3. Star formation models (fii = 0.1 and M tot = 10 10 M Q , unless noted otherwise) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 




Model 


(SFR) 


SFRi 


a 


[O/H]! 


iVtot 


L B 


M m / L B 


7 


Notes 






[Meyr- 1 ] 






[xlO 10 ] 


[xlO 9 L , B ] 


[Mq/L Qi b] 






Al 


Exp. decreasing 


0.93 


0.17 


0.18 


+0.3 


3.6 


3.6 


0.16 


-2.35 


(a) 


A2 




0.52 


0.09 


0.18 


-0.25 


2.1 


2.3 


1.56 


-2.35 


(a, b) 


A3 




0.69 


0.13 


0.18 


-0.9 


4.7 


1.6 


0.45 


-3 


(a) 


Bl 


Constant 


0.89 


0.89 


1.0 


+0.3 


3.4 


10 


0.03 


-2.35 




B2 




0.49 


0.49 


1.0 


-0.25 


2.1 


6.4 


0.58 


-2.35 


(b) 


B3 




0.68 


0.68 


1.0 


-0.9 


4.7 


3.7 


0.21 


-3 




C 


Exp. decreasing 


0.90 


0.07 


0.08 


+0.3 


11 


2.3 


0.32 


-2.35 


(c) 


D 


Lin. decreasing 


0.90 


0.60 


0.67 


+0.3 


3.5 


7.6 


0.10 


-2.35 




E 


Exp. increasing 


0.85 


3.06 


3.4 


+0.3 


3.3 


28 


0.03 


-2.35 




F 


Lin. increasing 


0.87 


1.74 


2.0 


+0.3 


3.3 


18 


0.04 


-2.35 





Notes: (a) r sfr = 5 Gyr. (b) (ii = 0.5. (c) r s f r = 4 Gyr 



tio of current and average past SFRs a, and the present- 
day mass averaged ISM oxygen abundance. The aver- 
age past SFR is roughly the same for all models ending 
at fjn = 0.1 and is (SFR) = 0.9 Moyr" 1 . In contrast, 
present-day SFRs range from SFRi = 0.07 to 3 M yr _1 
and in fact determine the contribution by young stars to 
the integrated light of the model galaxy. Current oxygen 
abundances predicted are [0/H]i <~ +0.25 and are mainly 
determined by /xi, the IMF, and the assumed mass limits 
for SNII (cf. Table 3). Column (7) gives the number of 
stars formed, Column (8) their total B band luminosity. 
Column (10) shows the resulting Hi mass-to-light ratio. 
Column (11) gives the slope of the IMF used, while Col- 
umn (11) refers to additional Notes. 



Let us consider the photometric evolution of the con- 
stant and exponentially decreasing SFR models in some 
more detail. Fig. 3 shows the evolution of the total num- 
ber of MS and post-MS stars for exponentially decreasing 
SFR model Al. The current total number of MS stars is 
roughly 4 • 10 10 , compared to ~ 10 8 post-MS stars (all 
phases) and ~ 3 • 10 4 for AGB stars only. The present-day 
mean luminosity of stars on the MS is ~ 10 _1 Lq com- 
pared to ~ 10 4 Lq for AGB stars. The current bolomet- 
ric galaxy luminosity is determined mainly by MS stars 
(Lms ~ 4 • 10 9 Lq). In particular, AGB stars (Lagb ~ 
10 8 Lq) are relatively unimportant. 

We show in Fig. 4 the B and /-band magnitudes of 
stars in different evolutionary phases (again for the expo- 
nential model). As for the bolometric galaxy luminosity, 



8 



L.B. van den Hoek et al.: Evolution of LSB galaxies 




5 10 15 5 10 15 



Age [Gyr] Age [Gyr] 

Fig. 3. Exponentially decreasing SFR model: total number of stars formed (left) and average stellar bolometric lu- 
minosity (right) vs galactic age for different evolutionary phases: MS (solid curve), RGB (dashed), HB (dot-dashed), 
EAGB (dotted), and AGB (dash-dotted). 




7 8 9 10 7 8 9 10 

Log Age [yr] Log Age [yr] 

Fig. 4. Photometric evolution of exponentially decreasing SFR (top panels) and constant SFR (bottom) models. Left 
panels: £>-mag contribution for stars in different evolutionary phases as in Fig. 3: MS (solid curve), RGB (dashed), 
HB (dot-dashed), EAGB (dotted), AGB (dash-dotted), and Total (thick solid). Right panels: Same as left panels but 
for /-mag 
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Fig. 5. Broadband colours vs galactic age for exponen- 
tially decreasing SFR (thick solid) and constant SFR 
(solid) models. Note the different magnitude scales in the 
upper two and bottom three panels, respectively. 

MS stars generally dominate in the -B-band. However, in 
the J-band, RGB and HB stars are nearly as important 
as MS stars, at least in later stages of galactic evolution. 
Due to the cooling of old, low-mass MS stars as well as 
the increasing contribution by RGB and HB stars with 
age, the current total /-band magnitude is considerably 
brighter than that in the B-band. This qualitative model 
behaviour is insensitive to the adopted star formation his- 
tory (for e- folding times larger than 3 — 5 Gyr) but instead 
is determined by the assumed IMF and the stellar input 
data used. Thus, constant star formation models exhibit a 
similar behaviour apart from being brighter by about one 
magnitude in all passbands at late times (cf. Fig. 4). 

Fig. 5 illustrates the sensitivity of broadband colours 
to the star formation history for constant and exponen- 
tially decreasing SFR models. The colours considered here 
redden with galactic age. In general, differences between 
colours such as for different SFR models are less than the 
variations of these colours with age for a given model. As- 
suming a galactic age of e.g. 8 instead of 14 Gyr has lim- 
ited effect on the resulting galaxy colours (e.g. less than 
0.1 mag in B — V), even though absolute magnitudes are 
substantially altered. Both age and extinction effects can 
result in substantial reddening of the colours of a stellar 
population in almost the same manner and it is difficult to 
disentangle their effects on the basis of photometry data 



alone. Clearly, galaxy colours alone are not well suited 
to discriminate between different SFR models, even when 
internal extinction is low and other reddening effects are 
negligible (see below). 



As our emphasis will be on a comparison between the 
evolution of HSB and LSB galaxies, we will first focus on 
the relatively well-known evolution of HSB galaxies, using 
a model chosen to emphasise the observed properties of 
HSB galaxies. Using this as a starting point, we explore 
other models to find the best model for LSB galaxies. 



Models using exponentially decreasing SFRs and differ- 
ent e-folding times for different Hubble types (e.g. Larson 
& Tinsley 1978; Guiderdoni & Rocca-Volmerange 1987; 
Kennicutt 1989; Bruzual & Chariot 1993; Fritze-v. Al- 
vensleben & Gerhard 1994) are usually found to be in 
good agreement with observations of HSB galaxies. There- 
fore we will mainly focus on the exponential SFR models 
discussed above. As the gas-richness is one of the obvious 
differences between HSB and LSB galaxies, we will first 
focus on models with \i\ =0.1 to describe HSB galaxies. 

Fig. 6 shows the results of the exponentially decreas- 
ing SFR model: SFR cx exp(— i/r s f r ) with r s f r = 5 Gyr 
ending at [i\ — 0.1 at t — t ev = 14 Gyr for several initial 
galaxy masses M g (t = 0) - 10 8 - 10 10 M@. The precise 
value of fj,\ was achieved by scaling the amplitude of the 
SFR accordingly. Galaxy colours and abundances are not 
affected by this scaling while absolute magnitudes and fi- 
nal gas masses scale with M g (t = 0) (as indicated in Fig. 
6). 

Present-day (B — V) and (R — /) colours for this model 
are 0.6 and 0.5, respectively. HSB galaxies with (B — V) <; 
0.6 and (R — I) 0.5 can be explained only when substan- 
tial amounts of internal dust extinction are incorporated. 
The reason is that even though a single stellar popula- 
tion born with metallicity Z — 0.02 may become as red as 
(B — V) ~ 1.1 and (R — I) <~ 0.8, the luminosity contribu- 
tion by young stellar populations (with ages less than a few 
Gyr) results in substantial blueing of the galaxy colours. 
For the same reason, values of (B — V) ^ 0.65 are not 
predicted by dust-free models independent of the adopted 
(exponential) SFR (as long as r ^ 4 Gyr) or IMF. There- 
fore, considerable reddening of the younger stellar popula- 
tions is required to explain the colours of HSB galaxies in 
this model. Figs. 6c and d show that both the [O/H] abun- 
dance ratios and present-day gas masses of HSB galaxies 
can be explained reasonably well for M g (t = 0) > 10 10 
Mq . These conclusions agree with results from many pre- 
vious models (with different stellar input data) for HSB 
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galaxies (e.g. Larson & Tinsley 1978; Guiderdoni & Rocca- 
Volmerange 1987). 

It is clear from Fig. 6 that this model is not able to 
explain the properties of the LSB galaxies in our sample. 
The model colours are too red while the predicted [O/H] 
abundances are too high. As LSB galaxies have higher 
gas-fractions we thus need to broaden the range of gas- 
fractions. 

5.2. LSB galaxies: \i\ ^ 0.1 models 

We now consider models with gas fractions \i\ = 0.025, 
0.1, 0.3, 0.5, 0.7, and 0.9. First we will examine the con- 
straints which individual properties such as abundance 
and gas-fraction impose on the models, before combining 
these in a final model. 

5.2.1. Colours and magnitudes 

Results for a model with initial mass M g (t = 0) = 10 10 
M Q are shown in Fig. 7. While the age distribution of the 
stellar populations is identical in all the models shown, 
present-day (B — V) and (R — I) colours are found to 
decrease by 0.2 and 0.1 mag, respectively, when going from 
models ending at Hi = 0.025 to 0.9. This bluing effect 
is due to the decrease of stellar metallicities for models 
ending at increasingly higher gas fractions. 

The LSB galaxies can be distinguished in two major 
groups by means of their colours. First, LSB galaxies with 
(B — V) ^ 0.5 which can best be modelled with expo- 
nentially decaying SFR models ending at values [i\ = 
0.5 — 0.7 at age 14 Gyr. These results can be shifted to- 
wards brighter or fainter B magnitudes by changing the 
initial gas masses accordingly. This leaves the resulting 
colours unaltered (cf. Fig. 6). 

Second, relatively blue LSB galaxies with (B — V) ^ 
0.4 cannot be fitted by exponentially decreasing SFR mod- 
els with r ^ 5 alone (assuming i cv = 14 Gyr), regardless 
of their current gas fraction \ii . There are several possible 
explanations: a relatively young "surplus" stellar popula- 
tion may influence the galaxy colours resulting from an 
underlying exponentially decreasing SFR. Alternatively, 
these LSB galaxies may have started forming stars re- 
cently (e.g. 5 — 8 Gyr instead of 14 Gyr, cf. Fig. 7a), or they 
may have more slowly decaying or constant SFRs which 
generally result in present-day colours of (B — V) ^ 0.4 
and (R — I) 0.35 (but the latter option appears to be 
excluded because of their abundances and gas fractions; 
see below). 

5.2.2. Abundances 

The observed range in [O /H] abundances for LSB galaxies 
is well explained by exponentially decaying SFR models 



with iii <; 0.3 (cf. Fig. 7c). Constant SFR models end- 
ing at /Ui <; 0.3 are also consistent as the abundances of 
elements predominantly produced in massive stars arc in 
general determined by the present-day gas fraction \ii and 
are insensitive to the detailed underlying star formation 
history (e.g. Tinsley 1980). However, these latter models 
appear to be excluded by the measured gas masses (see 
5.2.4 and 5.2.5) (for a closed box model). Though the 
abundances of metal-poor LSB galaxies with [O/H] 55 — 1 
can be well explained by models ending at /ii = 0.9, it is 
likely that the exponential models presented here are an 
over-simplification for these gas-dominated systems. It is 
more likely that they have experienced a low and sporadic 
star formation rate. 



5.2.3. Gas content 

Present-day gas masses observed in LSB galaxies can be 
well fitted by exponentially decaying SFR models ending 
at iii ~ 0.5 (Fig. 7d). This is consistent with the range 
of /ii k, 0.3 derived from the [O/H] data. Exponentially 
decreasing SFR models ending at /ii <; 0.6 are inconsistent 
with the observations unless we assume that LSB galaxies 
have started forming stars only a few Gyr ago. 

Fig. 8 displays the gas mass vs (B — V) for constant 
and exponentially decaying SFR models. Exponentially 
decreasing SFR models are able to explain simultaneously 
values of m ~ 0.5 ± 0.2, (B - V) £ 0.5, and M g k, 10 9 
Mq, as observed for the majority of the LSB galaxies in 
our sample. Over the entire range of possible gas fractions, 
constant (or increasing) SFR models are only able to 
explain the bluest galaxies. 

Fig. 8 demonstrates that exponentially decreasing SFR 
models ending at [ii <~ 0.5 are in good agreement with 
the observed hydrogen mass-to- light ratios Mm/ Lb, hi 
contrast to constant (or slowly decreasing) SFR models 
(expect for the bluest systems). 

The conclusion that constant and slowly decreasing 
SFR models are inconsistent with the observed photomet- 
ric and chemical properties of LSB galaxies is based on 
the assumption of negligible amounts of dust extinction in 
these systems. However, if a considerable fraction of the 
LSB spirals in our sample would suffer from extinctions of 
E(B-v) = 0.1 — 0.25, the Mm/ Lb ratios predicted by con- 
stant SFR models would increase by a factor 1.4 — 2.5 after 
correction for internal extinction. In this manner, constant 
and slowly decreasing SFR models ending at gas fractions 
Iii = 0.3 — 0.7 after ~ 14 Gyr could also explain relatively 
red LSB galaxies with (B — V) <; 0.5 and Mm/ Lb ratios 
as large as ~ 1.5. Even though we argued in Sect. 2 that 
internal extinction in LSB galaxies is unlikely to exceed 
E(b-v) ~ 0.1, knowing the dust content of LSB galax- 
ies is obviously of crucial importance in modelling their 
evolution. 
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Fig. 6. Exponentially decaying SFR model Al: photometric and chemical evolution results. Curves have been drawn 
for initial galaxy masses M s (t = 0) = 10 s , 10 9 , and 10 10 M . For the M g (t = 0) = 10 9 Af© model, arrows indicate 
evolution times of 1, 2, 4, 8, and 14 Gyr, respectively. Symbols refer to the following galaxy samples: face-on HSB 
spirals (dots; de Jong & van der Kruit 1994), LSB spiral galaxies (triangles; de Blok et al. 1995), and dwarf irregulars 
(crosses; Melisse and Israel 1994). Typical observational errors are shown in the bottom left of each panel. 



5.2.4. Summary: the best model 

To summarise, we find that exponentially decreasing SFR 
models ending at \i\ = 0.3 — 0.7 are in agreement with 
the colours, magnitudes, [O/H] abundances, gas contents, 
and mass-to-light ratios observed for LSB galaxies with 
(B - V) £ 0.45. Blue LSB galaxies with {B -V) <, 0.45 
cannot be fitted by exponentially decreasing SFR mod- 
els without an additional light contribution from a young 
stellar population. Alternatively, such LSB galaxies may 
have experienced constant SFRs, or may be much younger 
than 14 Gyr. 

5.3. HSB/LSB: Just an extinction effect? 

Since the models require an internal extinction in HSB 
galaxies up to E(b_v) ~ 0.5 (corresponding to Ay ~ 
1.5 mag), independent of the SFH, one could argue that 
extinction may be the main explanation for the difference 
in colour observed between LSB and HSB galaxies. 



As shown in this paper the models indicate that the un- 
derlying stellar populations in LSB and HSB galaxies are 
still distinctly different due to differences in the chemical 
evolution of these systems. Therefore, extinction effects, 
although important, cannot be the entire explanation for 
the colour differences observed. 

6. The effects of small star formation bursts 

6.1. Luminosity contribution by HII regions in LSB 
galaxies 

The presence of Hn regions in all the LSB galaxies in our 
sample suggests that recent star formation is a common 
phenomenon in these systems, and may thus influence the 
observed properties. In the previous section we found that 
for some LSB galaxies a simple exponential SFR was un- 
able to explain the observed properties and that an addi- 
tional "burst" of young stars was needed. Here we examine 
the effects of such bursts on the colours of LSB galaxies. 
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Fig. 7. Photometric and chemical evolution results in case of exponentially decaying SFR models ending at /ii = 
0.025, 0.3, 0.5, 0.7, and 0.9 for an initial galaxy mass of M g (t = 0) = 10 10 Mq (in Fig.7d we show also results for 
M g (t = 0) = 10 9 Mq). Arrows on top of the /ii = 0.025 model indicate evolution times of 1, 2, 4, 8, and 14 Gyr, 
respectively. Observational data as in Fig. 6 



To this end, we selected all Hn regions that could be iden- 
tified by eye, either from Ha or i?-band CCD images, and 
added their total luminosity in the /-band. We restrict 
ourselves to the /-band data as to be least affected by ex- 
tinction that may be present in or near Hn regions (Mc- 
Gaugh 1994). We define r\i as the ratio of this Hn region 
integrated luminosity and the total luminosity of the cor- 
responding LSB galaxy. As we are probably incomplete at 
low Hn region luminosities, this ratio gives a lower limit 
for the Hn region contribution to the total luminosity. 

In columns (1) and (2) of Table 4, we list the LSB 
galaxy identification and number of Hn regions selected. 
The number of Hn regions identified within individual 
LSB galaxy ranges from a few to ~ 25. For the ensem- 
ble of Hn regions in each LSB galaxy, we tabulate the 
absolute / magnitudes as well as the corresponding ratios 
rji of the Hn region integrated luminosity and total LSB 
galaxy luminosity, in columns (3) and (4). Mean (B — V) 
colours for the Hn regions and for the LSB galaxy as a 
whole are given in columns (5) and (6), respectively. 



For most of the LSB galaxies in our sample, the con- 
tribution of the Hn regions to the total light emitted by 
LSB galaxies does not exceed rjj = 0.05 — 0.1. However, 
as Hn regions are likely to contain an increased amount 
of dust (McGaugh 1994) the actual contributions may be 
higher by a factor of ~ 2 in /. Thus, the values of r/j pro- 
vide lower limits to the actual luminosity contributions of 
the Hn regions. For some LSB galaxies, e.g. F568-V1 and 
F577-V1, the Hn region contribution is found to be as high 
as rji = 0.2. These systems contain a only modest number 
of Hn regions so that their Hn regions on average may be 
larger and/or brighter than those present in several other 
LSB galaxies. 

Figure 10 shows the resulting Hn region contributions 
r\i for the SFR models discussed in Sect. 4. If we assume a 
maximum age run = 5 Myr for the Hn regions observed, 
then this implies that stars more massive than ~ 25 Mq 
are associated with these regions, according to the stellar 
evolution tracks from the Geneva group (see below). 

Corrections for dust extinction within the Hn regions 
and LSB galaxy as a whole, respectively, will shift the 
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Table 4. Observed broadband contributions by Hn regions in LSB 
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* values refer to R band magnitudes instead of /-band 

| uncertain due to contamination by fore- or background objects 
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observations in the directions as indicated in Fig. 10 (as- 
suming a mean Galactic extinction curve). From the / 
band observations we can conclude that the values of r]i 
predicted by smoothly varying SFR models are systemat- 
ically too low for most of the LSB galaxies in our sample. 

This is true in particular for F568-V1, F577-V1, and 
U628 for which values of rjj^ 0.2 — 0.25 suggest that star 
formation has been recently enhanced by factors ~ 5 — 10 
relative to the SFRs predicted by exponentially decreasing 
or constant SFR models. 



6. 2. Effects of small amplitude star formation bursts 

We investigate whether the observed values of i]i ~ 0.25 in 
the / band, as observed for several LSB galaxies discussed 
above, can be explained by small-amplitude bursts of star 
formation. 

We assume a Gaussian star formation burst profile 
with a given amplitude Af, and dispersion a. We follow 
its evolution during 1 Gyr with a time resolution of ~ 0.1 
Myr at time of burst maximum and of ~ 2 Myr at roughly 
5cr from burst maximum. We superimpose the star forma- 
tion burst on an exponential SFR model, as discussed in 
Sect. 5.1. 

We initially assume a galactic evolution time at burst 
maximum of % = 13 Gyr, burst maximum amplitude 
Ah = 8 Mq yr^ 1 , burst duration Aib = 5 Myr, max- 
imum Hn region lifetime thii = 5 Myr, and an initial 
galaxy mass of 10 10 M Q . For simplicity, we neglect any 
influence of the burst on the chemical evolution of the 
model galaxy. 



We show the result in Fig. 11a. At t = % = 13 Gyr 
the contribution by young stars increases rapidly. Simul- 
taneously, the (R — I) colours become significantly bluer. 
After burst maximum, the contribution by young stars 
decreases again and galaxy colours start to redden until 
the effect of the burst becomes negligible and colours and 
magnitudes evolve as prior to the burst. 

In this manner, a characteristic burst loop is completed 
as shown in Fig. 11a. The shape of this loop is determined 
by: i) the burst amplitude, ii) the extinction within the 
Hn regions, iii) the contribution by the old stellar popula- 
tion to the integrated galaxy light, iv) the duration (and 
profile) of the burst, v) the maximum time during which 
young stars produced by the Hn regions can still be dis- 
tinguished from the surrounding field population. 

Burst amplitude and IMF: Fig. 11a demonstrates that 
values of r/i ~ 0.25 can be well explained by bursts with 
amplitude A\, = 8 M Q yr _1 superimposed on an exponen- 
tially decreasing SFR model ending at /Ui = 0.1. Decreas- 
ing the burst amplitude by factors 2.5 and 5, respectively, 
results in the smaller loops shown in Fig. 11a. The ac- 
tual burst amplitude required to explain the observations 
depends on many quantities as described below 

Dust extinction: Fig. lib illustrates the effect of dust 
extinction on rjj and (R — I). A selective extinction of 
Eb-v = (0.25, 0.5, 1) mag for the Hn regions results in 
a reduction r\i by a factor (1.6, 2.5, 6.5) and a reddening 
Er_i = (0.14, 0.27, 0.56) mag, assuming a Galactic ex- 
tinction curve. For values of Eb-v ^ 0.5 mag, we find that 
the bluing effect on (R — I) by young massive stars formed 
during the burst is neutralised almost entirely by extinc- 
tion. For intense bursts reddening of the galaxy colours 



14 



L.B. van den Hoek et al.: Evolution of LSB galaxies 





Fig. 9. Evolution of Mm/ Lb vs (B — V) (left) and Mm/ Lb vs \i\ (right). Top panels: exponentially decreasing SFR 
models, bottom panels: constant SFR models. From bottom to top, curves are shown for models ending at increasing 
gas fractions /ii = 0.1, 0.3, 0.5, and 0.7, respectively. Arrows on top of the fix = 0.1 model indicate evolution times of 
1, 2, 4, 8, and 14 Gyr, respectively. Observational data as in Fig. 8 



may even occur. If variations in the mean extinction of 
the ensembles of Hn regions among LSB galaxies are small 
(e.g. less than a factor two), it is difficult to see how ex- 
tinction alone can provide an adequate explanation for the 
large variations observed in -qj. 

Global star formation history: From Fig. 11c it is clear 
that the burst effect on the galaxy magnitudes and colours 
increases when the contribution by the old stellar popu- 
lation is decreased. Thus, colours and magnitudes are less 
affected by bursts imposed on constant or even increas- 
ing SFRs compared to those imposed on exponentially 
decreasing SFR models: the smaller loop sizes just reflect 
that the mean age of the underlying stellar population is 
relatively low. 

Current gas fraction: Fig. lid demonstrates how the 
present-day gas fraction [i\ affects the effect of the burst. 



The luminosity contribution by the old stellar population 
decreases for increasing values of fii, and thus, as discussed 
above, we find that the effect of a star formation burst 
for galaxies with fj,i — 0.9 (i.e. unevolved systems) is as 
large as that of a ten times stronger burst for fix = 0.1 
(i.e. highly evolved systems). Thus, the burst amplitude 
required to explain the observations strongly depends on 
the present-day gas-to-total mass-ratio. 

Burst duration: Fig. lie shows that the duration of the 
burst affects the impact of the burst as well. For Atb = 1, 
5, and 10 Myr, the variation in rji is ~ 0.25 while the 
resulting galaxy (R — I) colours become bluer. Burst du- 
rations in excess of Atb ~ 5 Myr are unlikely since this 
would require dust extinctions Em—v) <^ 1 mag in order 
to provide agreement with the observed (R—I) colours (cf. 
Fig. lib). Such large extinction in Hn regions are probably 
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excluded by the observations (e.g. McGaugh 1994). Thus, 
relatively narrow burst profiles are needed to explain ex- 
treme values 77/ ~ 0.2 (as in F568-V1; cf. Table 4). 

Maximum age of H11 regions: Fig. 1 If shows that when 
the maximum lifetime of the H11 regions is increased from 
T~mi — 5 to 50 Myr, partial agreement with the obser- 
vations can be achieved without invoking star formation 
bursts. We note that the resulting H11 region contribu- 
tions do not increase linearly with thii as short lived mas- 
sive stars dominate the luminosity contribution of all stars 
formed during the past thii yr. 




7 

0.5 1 

B - V 




7 

0.5 1 

B — V 

Fig. 8. Comparison of predicted and observed gas masses 
vs (B — V). Top: exponentially decreasing SFR models. 
Bottom: constant SFR models. Models with SFR normal- 
isations according to m — 0.025, 0.1, 0.3, 0.5, 0.7, and 0.9 
are shown for an initial galaxy mass of M g (t = 0) = 10 10 
Mq. Arrows on top of the [i\ = 0.025 model indicate 
evolution times of 1, 2, 4, 8, and 14 Gyr, respectively. Ob- 
servational data as in Fig. 6 
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Fig. 10. Evolution of the Hn region integrated luminosity 
contribution for different SFR models ending at /ii = 0.1. 
Shown is the ratio of Hn region integrated and total lumi- 
nosity in the I band rjj vs (R — I). For the exponentially 
decreasing SFR model, solid triangles indicate evolution 
times of 1, 2, 4, 8, and 14 Gyr, respectively. Observational 
data on LSB galaxies is shown as open triangles. Typical 
errors in the data (upper right) and directions to which 
the data would shift after corrections for extinction (bot- 
tom left) are indicated. Top and bottom arrows indicate 
corrections for extinction within the Hn regions and the 
LSB galaxy for =0.1, respectively 

Values of thii ^ 50 Myr would mean that stars down to 
masses of m ion & 7 M Q would contribute to the Hn regions 
identified (e.g. Schaerer et al. 1993; Z — 0.001). Observa- 
tional estimates for m; on are usually in the range 10 — 15 
Mq (Wilcots 1994; Garcia- Vargas 1995) and correspond to 
thii i<5 15 Myr. Even though these values would imply that 
our adopted value of thii = 5 Myr is too low, this excludes 
extreme values of thii = 200 Myr which would be required 
to explain the observed range in r\i exclusively in terms 
of variations in thii and/or m\ m . Therefore, variations in 
thii (or equivalently mj on ) and/or extinction may provide 
an explanation only for part of the variations in the Hn re- 
gion integrated luminosity contributions observed among 
LSB galaxies. 

In conclusion, the observations suggest that small am- 
plitude, short bursts of star formation are important in 
at least several of the LSB galaxies for which accurate 
photometry data is available. Such recent episodes of en- 
hanced star formation may play an important role in af- 
fecting the colours of the blue LSB galaxies. 

6.3. Quantitative effect of small amplitude bursts on 
galaxy colours and magnitudes 

Table 5 lists the effect of a 5 Myr star formation burst 
on the galaxy integrated magnitudes and colours for var- 
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Table 5. Maximum changes induced on galaxy magnitudes and colours for M g = 10 10 M Q 
by a 5 Myr starburst 



Model 


A h 


AB 


AI 


A(B — V) 


A(R-I) 




Notes 




[Mq yr" 1 ] 


[mag] 


[mag] 


[mag] 


[mag] 






Al+burst 


0.8 


-0.34 


-0.14 


-0.18 


-0.04 


SFRi 


= 0.17 Moyr" 1 


» 


1.6 


-0.58 


-0.30 


-0.28 


-0.10 






V 


3.0 


-1.05 


-0.42 


-0.43 


-0.17 




* 


•)■) 


8.0 


-1.72 


-0.53 


-0.56 


-0.26 




* 


Bl+burst 


0.8 


-0.14 


-0.06 


-0.08 


-0.02 


SFRi 


= 0.89 M yr _1 




1.6 


-0.32 


-0.21 


-0.11 


-0.07 








3.0 


-0.56 


-0.32 


-0.21 


-0.11 




* 


!J 


8.0 


-0.94 


-0.41 


-0.33 


-0.14 




* 



* E( B _v) 2; 0.3 mag in Hn regions is required to provide agreement with observations 



Table 6. Model SFR 
and gas fractions 



(1) 


(2) 


Ml 


A 




(Meyr- 1 ) 


0.1 


0.18 


0.3 


0.13 


0.5 


0.09 


0.7 


0.06 


0.9 


0.02 



ious burst amplitudes. For an exponentially decreasing 
SFR (model Al in Table 3; assuming M g (0) = 10 10 M , 
//i =0.1, Salpeter IMF), we find that a burst with ampli- 
tude A^ = 0.8 Mq yr -1 results in maximum colour vari- 
ations A(B - V) and A(R - I) of -0.18 and -0.04 mag, 
respectively. The effect of increasing the burst amplitude 
by a factor 10 to A^ = 8 M Q yr" 1 results in corresponding 
shifts of —0.56 and —0.26 mag, respectively. This effect is 
similar to that when the initial galaxy mass is reduced by 
a factor ten while leaving the burst amplitude unaltered 
(i.e. for. M g (0) = 10 9 M and A h = 0.8 M yr" 1 ). 

For bursts superimposed on exponentially decreasing 
SFRs models ending at //i = 0.1, the colour and magni- 
tude shifts predicted are consistent with the observations 
in case of burst amplitudes A D ^ 3 M yr" 1 , assuming 
a typical extinction of E/ B _ v \ = 0.3 mag in the Hn re- 
gions. In fact, the effect of the burst is determined mainly 
by the total luminosity of the young stellar populations 
formed according to the continuous SFR during the last 
Gyr or so. Since the amplitude of the SFR scales with 
(1 — /ti) for models ending at different gas fractions /ti 
the impact of the burst for models ending at //i = 0.5 is 
about twice that given in Table 5. Similarly, for models 
ending at //i <; 0.9, the burst effect becomes roughly ten 
times stronger compared to the /ii = 0.1 case (cf. Fig. 11). 
The effect of the burst is substantially reduced when go- 
ing from exponentially decreasing to constant SFRs (cf. 
Table 5). 



7. Present-day star formation rates in LSB 
galaxies 

7.1. Theoretical star formation rates 

7.1.1. Continuous SFR 

Theoretical SFRs can be derived from the models dis- 
cussed in the previous section (exponentially decreasing 
SFR, Salpeter IMF), using: 
M tot 

SFR-^AG^)— Moyr" 1 (3) 

where A(fj,i) is the model SFR amplitude required to end 
at a gas fraction //i at a galactic evolution time of 14 Gyr 
(assuming an initial mass of 10 10 M ) and M tot the total 
galaxy mass as obtained from Mm and H\. 

We list the required values for A to have the models 
end at certain values of //i in Table 6. Using these values 
we find that LSB galaxies show present-day SFRs (without 
bursts) between 0.01 <> SFR cont <> 0.15M Q yr- 1 . For a 
typical LSB galaxy (i.e. M tot = 10 10 M Q , /*i = 0.5) we 
estimate SFR cont = 0.1 Moyr" 1 . 

7.1.2. Burst SFR 

As we discussed in Sect. 6, the effect of a starburst on 
the galaxy integrated magnitudes and colours depends on 
many quantities. However, assuming a modest Hn region 
extinction of E(b-v) — 0.25 mag and lifetime ran = 20 
Myr, crude estimates for the maximum burst amplitude 
can be derived from: 

SFR^« 20*^(^-1) MQyr" 1 (4) 

We find that /-band contributions of tji = 0.2 are best 
explained by a 5 Myr burst with amplitude SFR burst = 0.8 
M Q yr _1 , assuming present-day values of Mm = 2-10 9 M Q 
and // rot = 0.5. 

Theoretical estimates for the total current SFRs in 
LSB galaxies are found from SFR tot = SFR cont + SFR burst 
and are listed in the last two columns of Table 7. Present- 
day SFRs for LSB galaxies experiencing small amplitude 
bursts range from SFR tot - 0.02 to 0.8 M yr" 1 . 
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Table 7. Observational estimates of present-day SFRs [M Q yr 1 ] 


in LSB 


galaxies* 








(!) 


(2) 


(3) 


/ A \ 

( 4 ) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Name 


1J1&L. 


thj 


R 27 


P rr r 


Hi 


IVlHI 


(am) 


SFR? X 


SFRi 


g FR cont 


SFR tot 




[Mpc] 


[mag] 


[kpc] 


[kpc] 


[mag arcsec -2 ] 


[Mq] 


[M e p C - 




[Moyr- 1 ] 




F561-1 


47 


14.7 


6.6 


7.6 


23.2 


8.91 


4.5 


0.05 


0.09 


0.013 


0.081 


F563-1 


34 


16.1 


10.2 


16.3 


26.4 


9.19 


1.0 


(0.03) 


(0.02) 


0.16 


1.7f 


F563-V1 


38 


15.8 


4.6 


4.8 


24.0 


8.48 


4.2 


0.02 


0.03 


0.005 


0.024 


F567-2 


56 


15.9 


8.4 


10.6 


24.6 


9.09 


3.5 


0.04 


0.06 


0.029 


0.15 


F568-1 


64 


15.0 


9.6 


11.5 


23.7 


9.35 


5.4 


0.08 


0.18 


0.12 


0.31 


F568-3 


58 


14.7 


8.7 


11.4 


23.3 


9.20 


3.9 


0.09 


0.14 


0.06 


0.33 


F568-V1 


60 


15.7 


8.4 


10.7 


24.3 


9.14 


3.8 


0.05 


0.07 


0.05 


0.58 


U128 


48 


13.5 


18.2 


21.4 


24.1 


9.55 


(2.0) 


(0.16) 


(0.21) 


0.14 


0.82 


U1230 


40 


14.3 


12.0 


18.8 


24.5 


9.51 


(4.3) 


(0.11) 


(0.15) 


0.05 


0.39 



Notes: * theoretical values in columns (11) and (12) repeat those in columns (9) and (11) from Table 4; f probably contaminated 
by field galaxy. Values between parentheses are uncertain. 



7. 2. Empirical star formation rates 

We have also derived current SFRs for LSB galaxies us- 
ing the empirical method presented by Ryder & Dopita 
(1994; hereafter RD) based on CCD surface photometry 
of galactic disks. These authors found a relationship be- 
tween the local Ha and /-band surface brightness in the 
disks of a sample of 34 of southern spiral galaxies. From 
this relation, RD derived a constraint on the present-day 
SFR integrated over the entire stellar mass spectrum as: 

log SFRi = -0.26/2/ + 0.92 log cthi + 5.3 

Mopc^Gyr" 1 (5) 

where \ii is the /-band surface brightness and cm is the 
global mean Hi surface density [M©pc -2 ] within the star- 
forming disk. The relation between SFRi and fii is nor- 
malised by a term related to the mean surface density ohi 
and by a constant which is partly related to the conversion 
of the massive star formation rate to total SFR depend- 
ing on the adopted IMF (cf. Kennicutt 1983). It is unclear 
from the RD sample whether the relation is valid also for 
the lowest (stellar) /-band surface densities that are ob- 
served among LSB galaxies. However, since this relation 
holds over a wide range in surface brightness and massive 
star formation in the disks of spirals it appears rather in- 
sensitive to galactic dynamics, extinction, and molecular 
gas content (RD), we expect this relation to be valid also 
in case of LSB galaxies. At the faintest surface bright- 
nesses (i.e. fii ^ 25.6 mag arcsec -2 ), the relation may be 
flattening off although the effects of sky subtraction and 
small number statistics leave this open to question. If flat- 
tening indeed occurs, Eq. (5) provides lower limits to the 
actual SFRs in LSB galaxies. 

Using Eq. (3), we estimate global present-day SFRs 
for all LSB galaxies in our sample with measured /-band 
magnitudes and related data. For these LSB galaxies we 
list the distance, apparent /-band magnitude, and radius 
/?27 of the 27 mag arcsec -2 £>-band isophote, in columns 
(2) to (4) in Table 7. This radius corresponds to the optical 



edge of the LSB galaxy and is more representative for the 
radius within which the old disk stellar population in LSB 
galaxies is contained than is R25 as used by RD for HSB 
galaxies. Accordingly, we define an effective I band surface 
brightness as: 

nf = m/ + 2.51og(7r/? 2 7 ) (6) 

and use this in Eq. (5) . We tabulate the outermost radius 
of the measured Hi rotation curve /?hi, effective / band 
surface brightness /if*, total Hi mass derived within Rui, 
and the mean global surface Hi density cthi in columns (5) 
to (8), respectively. For a consistent comparison between 
effective surface brightness and global Hi surface density, 
one ought to measure them out to the same radius (e.g. 
/?27)- However, since Mm has been measured within Rm, 
we expect the former values to be more representative of 
the average Hi surface density in the star forming part of 
the disk. The mean global Hi surface densities cthi derived 
using /?hi vary between 2 and 5.5 M Q pc -2 and are sub- 
stantially smaller (i.e. by 20 — 60 %) than those derived 
using JS27 instead. 

The empirically derived mean present-day SFRs for 
the LSB galaxies from our sample range from about SFRi = 
0.02 to - 0.2 M yr -1 (cf. column 10 of Table 7). Er- 
rors arising from the Hi normalisation are estimated to be 
within a factor of two. This is illustrated when the same 
SFRs are derived assuming a fixed Hi surface density of 2 
M p C - 2 for all LSB galaxies (cf. SFRf x in Table 7). 

The empirically derived current SFRs in LSB galax- 
ies are in good agreement with the theoretically derived 
SFRs ranging from SFR cont - 0.01 to 0.15 A^yr -1 (cf. 
Table 7). As discussed before, the theoretically derived 
present-day SFRs of individual LSB galaxies lie proba- 
bly between SFR cont and SFR tot where the latter val- 
ues include the contribution of small amplitude bursts 
(SFR tot = SFR cont +SFR burst ). In several cases, the val- 
ues of SFR tot are considerably larger than the empirical 
values which suggests that the burst contribution is over- 
estimated by the models and/or that the SFRs derived 



18 



L.B. van den Hoek et al.: Evolution of LSB galaxies 



0.2 



0,1 - 



. \ a) amplitude . 

■ VA1V - 

i , i , i , i 


\ 

i , i 


\ b) extinction 

\ .". A 

\ ! 

1 , 1 




ii 1 i 


1 1 1 

d) gas fraction 






A 


. ^^^^^ . 




\ A 














i , i , i , i 


i , i 




,1,1,1 


\ 1 ■. ■• ■ 

. \ ; '. QHIIage 

: \\\\ A 1 
- \ "A 

: \ ".^ 



0.2 





0.3 



_1 ; b) effect of mean extinction within Hn regions: 
0.25 (solid curve), 0.5 (dotted), and 1 (dashed); 



o - 
i i 

0,2 0,4 0,6 0.2 0.4 0.6 
R-I 

Fig. 11. Impact of small amplitude bursts on the evolu- 
tion of the Hn region integrated luminosity contribution 
in the 7-band r/i vs (R — I). Unless stated otherwise, we 
plot the exponentially decreasing SFR + burst model as- 
suming A h = 8 M & yr -1 , Ai b = 5 Myr, r H n = 5 Myr, 
and Eb-v = mag. a) effect of varying burst ampli- 
tude: Ab = 8 (solid curve), 4 (dotted), and 1.6 (dashed) 
M yr- 
Eb-v = 

c) effect of global star formation history: exponentially de- 
creasing SFR (solid curve), constant (dotted), and linearly 
increasing (dashed); d) effect of actual gas-to-total mass 
ratio: = 0.1 (dashed curve), 0.5 (solid), and 0.7 (dot- 
ted); e) effect of burst duration: Aib = 1 (dashed curve), 
5 (solid), and 10 (dotted) Myr; f) effect of Hn region age: 
thii = 2 (solid curve), 10 (dotted), and 50 (dashed; no 
burst shown) Myr. Symbols have the same meaning as in 
Fig. 10. 



empirically do not trace well local enhancements of star 
formation at the faint surface brightnesses of LSB galax- 
ies. 

The present-day SFRs in LSB galaxies are thus found 
to be considerably below the ^5 — 10 M Q yr _1 derived 
for their HSB counterparts (e.g. Kennicutt 1992 and ref- 
erences therein) but significantly larger than the <~ 0.001 
M yr _1 observed typically in dwarf irregular galaxies (e.g. 
Hunter & Gallagher 1986). 



8. Discussion 

In the previous sections we have focussed on the difference 
between HSB galaxies and the type of LSB galaxies we 
have in our sample. As noted in the introduction (and 
also shown in Fig. 1) HSB and LSB galaxies, however, do 
not form separate groups. In fact there is a broad range 
of galaxy properties and the purpose of the discussion is 
to clarify what role the star formation history plays in 
determining the final character of a galaxy. 

In general, there appears a trend along the Hubble se- 
quence from rapidly decaying SFRs for early type galax- 
ies to constant or even increasing SFRs for dwarf irregu- 
lar galaxies (see e.g. reviews by Sandage 1986; Kennicutt 
1992). On average, the observed trend corresponds to a de- 
crease of the ratio of mean past to present SFR along the 
Hubble sequence. Most of the LSB galaxies in our sample 
belong to the group of late-type galaxies for which expo- 
nentially decreasing SFR models are in good agreement 
with the observations. The remaining LSB galaxies, for 
which slowly decreasing or constant (sporadic) SFR mod- 
els are more appropriate, belong to a group of galaxies 
with properties intermediate to those of disk galaxies with 
weak or absent spiral arms and Sm/Im galaxies. Thus, in 
general LSB galaxies comply well with the observed trend 
of SFR variation with Hubble type. 

The "burst" scenario discussed in section 6 indicates 
that current star formation in virtually all the LSB galax- 
ies in our sample is local both in time and space and sug- 
gests that sporadic star formation has been a continuous 
process from the time star formation started in the disks 
of LSB galaxies. 

The low star formation rates of ~0.1 M Q yr _1 derived 
for LSB galaxies as well as the local nature of the star 
formation in these systems are consistent with the idea of 
a critical threshold for the onset of global star formation in 
disk galaxies (e.g. Skillman 1987; Kennicutt 1989; Davies 
1990). 

Even though LSB galaxies contain large amounts of 
gas, only very limited amounts participate in the process 
of star formation. If we assume that LSB galaxies main- 
tain their current star formation rate of ~0.1 M Q yr _1 , 
their typical present-day amount of gas ~ M g = 2.5 10 9 
Mq will be consumed within r gas ^30 Gyr (for a recycled 
fraction of 25%). Such gas consumption times for LSB 
galaxies are much larger than a Hubble time (e.g. Ro- 
manishin 1980). For comparison, T gas ^2 — 4 Gyr in HSB 
galaxies, assuming typically M g ~ 10 10 M Q and SFRi ~ 5 
M Q yr _1 , which implies that HSB galaxies will run out of 
gas soon (see Kennicutt 1992). 

The presence of an old stellar population in many late- 
type LSB galaxies as indicated by their optical colours 
(e.g. vdH93; dB95) and as confirmed by the results from 
the modelling suggests that LSB galaxies roughly follow 
the same evolutionary history as HSB galaxies, but at a 
much lower rate. 
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The models suggest that very slowly decreasing star 
formation rates with e-folding times much larger than ~ 
5 Gyr probably can be ruled out for the majority of the 
LSB galaxies examined here. The observed colours are too 
blue, provided extinction is indeed negligible. This result is 
insensitive to the possible occurrence of infall or accretion 
of gas. Similarly the models indicate that the dominant 
stellar population in galaxies with (B — V) > 0.5 mag 
cannot be as young as 10 or 5 Gyr (cf. Fig. 5). In contrast 
the models can not rule out a dominant luminosity contri- 
bution by stellar populations significantly younger than 5 
- 10 Gyr for LSB galaxies with (B - V) w 0.4. 

First this indicates that the mean age of the stellar 
population in most LSB and HSB galaxies is similar even 
though the disks of LSB galaxies are in a relative early 
evolutionary stage. Although we have not explored an en- 
tire range of star formation e-folding times, values much 
below or in excess of ~ 5 Gyr are not very likely. Smaller 
values would increase the colour contribution from the un- 
derlying older stellar population and hence produce red- 
der colours than observed, while much larger values cor- 
respond to slowly decreasing or almost constant star for- 
mation rates implying a relatively larger contribution of 
the younger stellar population. 

Secondly, the combined effect of extinction and metal- 
licity on galaxy colours is sufficient to explain the colour 
differences observed between LSB galaxies and HSB galax- 
ies. Since the amount of extinction depends strongly on 
the dust content, which in turn is coupled to the heavy 
element abundances in the ISM (see Sect. 2), metallic- 
ity probably is the main quantity determining the colour 
differences observed between LSB and HSB galaxies. In 
this manner the much lower rate of star formation in LSB 
galaxies, implying lower metallicities and dust contents, 
indirectly determines the blue colours of LSB galaxies 
compared to HSB galaxies. 

9. Summary 

We have examined the star formation histories of LSB 
galaxies using models which take into account both the 
photometric and the chemical evolution of the galaxies. 
For the majority of the LSB galaxies in our sample, ob- 
served UBVRI magnitudes, [O/H] abundances, gas masses 
and fractions, and Hi mass-to-light ratios are best ex- 
plained by an exponentially decreasing global SFR ending 
at a present-day gas-to-total mass-ratio of [i\ ~ 0.5. When 
small amplitude bursts are involved to decrease the pre- 
dicted M g /L ratios, models ending at « 0.7 may also 
apply. In addition to exponentially decreasing SFR mod- 
els, ~ 15% of the LSB galaxies require modest amounts 
of internal extinction E(b-v) ^0.1 mag to explain the 
relatively red colours of (B — V) ~ 0.6 of these systems. 

A substantial fraction (~ 35%) of the LSB galaxies 
in our sample have colours (B — V) ^ 0.45 mag and 



exhibit properties similar to those resulting from expo- 
nentially decreasing SFR models at evolution times of 
^5—10 Gyr. Alternatively, recent episodes of enhanced 
star formation superimposed on exponentially decreasing 
SFR models may provide an adequate explanation for the 
colours of these systems (see Sect. 6). Recent star forma- 
tion is observed, at least at low levels, in essentially all 
the LSB galaxies in our sample. Hence, it seems justi- 
fied to assume that the disks in LSB galaxies experienced 
continuous (i.e. frequent small amplitude bursts of) star 
formation, at least during the last few Gyr. 

There is nothing peculiar about the evolution of LSB 
galaxies. Broadly speaking their evolution proceeds like 
that of HSB galaxies, but at a much lower rate. 
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